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SUMMARY 


This is the Final Report for Contract NAG5-1925, that consisted of experiment 
design, for possible use by the space science mission called WISP ( Waves in 
Space Plasma) • This mission is under study by the Canadian Space Agency and by 
NASA. Two WISP configurations are contemplated, under the name of BICEPS: one is 
called BOLAS, and the other WISPRS. Both these configurations are meant to per- 
form bistatic sounding of the ionosphere, at a height close to ( about 350 

Km ) f with a pair of satellites, either tethered or in free flight. 

INVESTIGATION A ( with Y.L. Alpert as P.I.) addressed the subject of parametric 
decay effects, expected to arise in a magnetoplasma under the influence of high- 
intensity HF fields. Criteria were formulated that could be used in searching for 
parametric instabilities and of electric fields in the ionosphere and magnetosphere, 
by in-situ satellites, such as the BICEPS pair. 

INVESTIGATION B ( with M.D.Grossi as P.I.) addressed the bistatic measurement^ by 
the BICEPS pair, of ionospheric features, such as large-scale and small-scale dis- 
turbances, Travelling Ionospheric Disturbances, electron density irregularities, 
spread-F phenomena, etc. These measurements by BICEPS could be correlated with the 
waveform distortion and degradation experienced by microwave links from geosynchro- 
nous height to ground, such as the ACTS satellite, expected to radiate pulses as 
short as 1 nanosecond in the band 20 to 30 GHz. These links are transionospheric and 
propagate e.m. waves in the volume of the ionosphere where BICEPS operates. It will 
be possible, therefore, to correlate the two classes of measurements, and learn the 
causative mechanisms that are resposible for the time-spread and frequency-spread 
nature of communications waveforms at microwave, in geosynchronous height to ground 
paths . 
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Abstract 


The essential results of this study, namely of the analytical and 
mainly of the detailed numerical calculations of some parametric de- 
cay effects, arising in a magnetoplasma under the influence of a HF, 
sufficiently strong electric field E = E p -cosut, are given in this final re- 
port. The resonance branches and the VLF parametric resonances are 
calculated in the ionosphere at altitudes Z = 200, 300400 km. The 
values of E p in the resonance regions were evaluated in the kinetic 
approximation. It is shown that the angle dependence u>i(0, E p / 0) 
of the ELF (0 < u) < ft#) resonance branch is close to the cosine law. 

This is in contrast with the earlier published results and with the angle 
dependence u>i(0, E = 0). This effect and the other discussed here 
dependencies may be used for search of parametric instabilities and of 
the electric field in the ionosphere and magnetosphere, especially by 
experiments in situ on satellites. 

I. Introduction 

The creation of the non-linear theory of parametric decay effects , pro- 
duced in a magnetoplasma under the action of a high frequency electric field 
E p cos u)£t } where loe is much larger or comparable with the Langmuir plas- 
ma frequency u> 0 , began about 30 years ago. See, for example, the papers 
by Silin and Aliev & Silin, 1964, 1965, 1966, [l to 3] ; Du-Bois Goldman, 
1965, 1967 [4, 5], Nishikawa, 1967, 1968 [6, 7], and the comprehensive mono- 
graph of Silin, 1973, [8]. However, until now these phenomena were a very 
little studied in the free plasma surrounding the Earth, in the ionosphere and 
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magnetosphere. Some causes of this situation may be the following. 

Such kind of experiments done in situ by satellites should be accompanied 
by delicate manipulations between the source of the electromagnetic field (the 
transmitter) and the receiver. The expected spectrum of the excited oscilla- 
tions is rather various. Thus, the experimental research and the discovery of 
these phenomena is not a simple task. On the other hand, the general theory 
predicts various non-linear parametric phenomena, which should be created 
in the ionosphere and magnetosphere. However, to find them, particularly 
on satellites and for successful theoretical treatment of experimental data, 
the knowledge of the precise conditions of the experiment, of the basic char- 
acteristics of the magnetoplasma, of the values of E p and u>e, etc. is required. 

To perform such experiments and to discover the parametric phenomena 
in the ionosphere and magnetosphere is important. It will be an essential 
contribution to the knowledge of the ionosphere and to the plasma physics in 
general. The parametric effects will be separated, divorced from many other 
linear and non-linear phenomena because of their special character. This will 
also help to understand other peculiar effects recorded in the ionosphere and 
magnetosphere which have not been explained until now. On the other hand, 
the theoretical treatment of these experiments will promote the development 
of the non-linear parametric theory of a magnetoplasma. Namely the noted 
reasons and the planned by NASA proper missions induced the theoretical 
calculations presented here. They are based on the general theory developed 
earlier (see [8]). 
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II. Some initial data and characteristics of the plasma 


The oscillations of the electrons of a plasma under the action of an alter- 
nating sufficiently strong electric field E p cosojEt are changing essentially the 
electromagnetic oscillatory behavior of the media. This happens when the 
velocity of the oscillation of the electrons te || $ becomes much larger or, 
in some cases, only comparable with the thermal velocity v e of the electrons 
when the frequency u>e is comparable with the electron Langmuir frequency 
u? 0 = (47rjV • e 2 /m) 1 ^ 2 . This condition which has the form 


e E p 

r E = >■ ve 

m u>e 


- ( 2kT '^ 1/2 


V m 




( 1 ) 


brings to the basic inequality of the phenomena discussed here 

El 

-£■ > NkT , 

4x 


( 2 ) 


where T is the temperature of the plasma particles. It means that the energy 
density of the electric field should be much larger than the thermal energy 
density of the electrons, i.e. larger of their pressure. In such conditions, 
a new kind of instabilities and resonances, and other peculiar properties of 
the spectra of the excited waves and oscillations appear in the plasma. This 
mechanism is the so called parametric decay instability and resonance of the 
plasma in the strong pump waves E = E p cos ust. The necessity of strong 
waves, i.e. of large amplitudes E v of the electric field, and also of large 
frequencies of the pumped wave u E u> 0 , is directly seen from (1) and 

(2). However, it is shown below, that in some cases, namely in the resonance 
region a >e ~ u>o, the amplitude of the electric field may be even much smaller 
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than, let us say, the characteristic amplitude Epo of the parametric non-linear 
effect which may be determined as 

E p0 = {4 tNkT') 1/2 . (3) 

Another important characteristic value which appears in the theory of the 
parametric non- linear effects in a plasma is the following: 

p p = k ■ r E = — —j = 1.91 • 10 7 • • u ■ n(u>). (4) 

A m u> E 

In (4) k — 2tt/A = w • n/c is the wave number, c is the velocity of e.m. waves 
in free space, n = n(u>) is the index of refraction of the local region of the 
plasma, and A is the wave length of the waves. 

In this study we are mainly interested in the very low - VLF (fi# < u> < 

fij,) and low - LF (fl/, < u < Ub) frequency bands. But some data are also 

given for the extra low ELF (0 < u> < fig) and high frequency HF (u>b < 

u) — ¥ uq...) bands. 1 , where Mb = (e • B)/(c- M e jj) and ub — (e- B)/(c-m) are 

— # 

the ion and electron angular gyro- frequencies and B is the constant magnetic 
field of the Earth. The ion Langmuir frequency u> ~ fl 0 = {4irN ■e 2 /M e ff) 1 ^ 2 
( M e fj is the effective mass of the ions) belongs to the LF frequency band. 
The results of calculations presented here were done for the ionosphere at 
altitudes Z = 200, 300, and 400 km. This region is described by the values 
of the parameters given in the following table, where besides the ones used 
above, some other characteristics of the plasma can be found. These values 
were used in the numerical calculations given in the next section. 

1 The classification of the frequency bands used here was introduced in [9]. 
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Characteristics of the Ionosphere 


Z, km 

N, cm~ 3 

"e., s _1 

ESI 

v e , cm/s 

Pi 

v„ cm/s 

w 0 , Hz 

200 

300 

400 

5.0 • 10 5 
1.8- 10 6 
1.5- 10 6 

2.0 • 10 3 
3.0 • 10 3 
1.5- 10 3 

1000 

1500 

1800 

1.1 • 10 7 
2.0 • 10 7 
2.8 • 10 7 

12 

5 

3 

7.4 ■ 10 4 
2.1 • 10 5 
3.8 • 10 5 

4.0 • 10 7 
7.6 • 10 7 
6.9 • 10 7 

Z, km 

f2o, Hz 

n L , Hz 

n B , Hz 



V A , m/s 

p c , cm 


200 

300 

400 

2.7 • 10 s 
7.9 • 10 s 
9.3 • 10 s 


5.3 • 10 4 
7.9- 10 4 
9.8 ■ 10 4 

364 

834 

1327 

■ 1 S ■ 

i ^ 



1.35 

2.60 

3.83 

Z, km 

Pi, cm 

4 nNnT e , 
erg /cm 3 

Ep o» 

CGSE 

E p o, 
V/m 

D e , cm 

D\, cm 2 



200 

300 

400 

203 

252 

286 

8.67 • 10~ 7 

4.68 • lO" 6 
4.68 • 10~ 6 

9.31 • 10~ 4 
2.16 - 10" 3 
2.16 ■ 10" 3 

27.9 

64.9 
64.9 

0.198 

0.186 

0.286 

3.78 • 10~ 2 
3.46 • lO" 2 
8.21 ■ lO" 2 



In addition to the notations given above, in this table Epo is the charac- 
teristic amplitude of the electric field (see (3)), D t is the Debye length of the 
electrons, D 3 = nT/AnNe 3 , the ions effective mass 


\/f Jl/f Mfj t 


(5) 


Hu, Mu are, respectively, the densities and masses of the ion constituents of 
the plasma, v, = yj2nT / M e jj is the thermal velocity of the ions, the subscript 
Hi is the notation of the protons, u>i is the low-hybrid frequency, p, = Vi/Qg, 
Pe/wg are the ion and electron Larmor radii, fi 0 = (4x7Vf/Af e //) 1 ^ 2 is the 
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ion Lengmuir frequency, fls and u >b are the circular gyro- frequencies of the 
ions and electrons, V4 = cjn 4 and n a = flo/Dg are the velocity and index 
of refraction of the Alfven waves. 

We consider here a quasi-neutral N e = TV, and isothermal T e = 7) plasma, 
i.e. that the energy densities and the Debye lengths of the electrons and ions 
are equal, 

N e tiT e ~ NiitTi, kD e ~ kDj. _ (6) 


III. Theoretical equations, formulas, results of numerical calculations 

The following theoretical results in this section give a qualitative and an 
approximate quantitatively representation of the behavior of the e.m. plas- 
ma oscillations and waves, excited parametrically in the ionosphere under 
the action of the electric field E p • cos ugi. 

The equations used here are based on different asymptotic solutions of 
the dispersion equation. For example, to calculate the resonance branches, 
the cold plasma approximation was mostly used. I.e. the values «,• = k • v, 
and u>d = k • D e , which appear in the kinetic approximation, due to the 
influence of the temperature of the particles, are absent in the formulas. 
Nevertheless, in the ELF and VLF frequency bands these values become 
comparable with the frequencies of the parametric waves and they should, 
and were, taken into account (see below). The cold plasma approximation 
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can be also used by calculation of the frequencies and growth rates of the 
waves in the VLF resonance region. However, the amplitudes E p of the elec- 
tric field can be, and were, evaluated by the kinetic approximation. In some 
cases, the isotropic plasma approximation was used. I.e. the influence of the 
constant magnetic field was not completely taken into account, and this may 
be warranted since wo u>b- 


III-l. Resonance branches of the ionosphere 


The dispersion equation of a magnetoplasma, not taking into account 
the thermal velocities of the electrons and ions, i.e. in the cold plasma 
approximation, can be conveniently represented as the following: 

D(u>, 0) = oj 4 (u>q — w 2 )(fi£ — <x 2 ) + u> 4 ($Iq(ub — u; 2 ) + u>q(Q.^ — a; 2 )) sin 2 0 — 
— a> 2 (u; o + — u; 2 )(f2^ — u; 2 ) cos 2 0 + (7) 

-I- (l — JliPv)) ' (“B cos2 ® — w2 )(^b cos2 ® ~ w 2 ) = 0, 


when u)£ > w 0 , U E > wb, and 
2irpi 2 5r/) e 

Ax ’ A ± < 


V i _ 

&B 



— < 

1 - s 

, < 

1 — s — 

v *\\ 

UJ 

V *\\ 

UJ 


( 8 ) 


where s = 0, ±1, ±2, .., p e , p,, and p p were defined above, 


v n = 


2 irc u*x 
Aj_ = = Z7T . 


UJ 


( 9 ) 


7l|| 71 

and v*||, u*x are the longitudinal and transverse phase velocities of the waves 
with frequency u>, and ny and nx are their longitudinal and transverse in- 
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dexes of refraction (see [8]). 


The properties of the roots of the equation D( u>, ©) = 0 are directly vis- 
ible. It estimates four wave branches u> 4 (0), s = 1 to 4. Only one of them, 
let us denote it u> 2 (©), is of pure electric field origin (see (11) and (12)). The 
influence of the electric field on the other resonance branches becomes clear 
only by analyzing the results of numerical calculations of the angle depen- 
dencies of u; 4 (0), where © is the angle between the wave vector k and the 
constant magnetic field B of the Earth. 

When 0 = 0, the frequencies a; 4 (0) are the following: 

W 1 = w 3 = ( 10 ) 

w 2,4 = \ ((Wo + flo) T y/H + n o) 2 - > J* ~ ^2 _ ^2 ? 

where 

tij, = (l — ^o(Pp)) > Pp = ~~f' ' ~r- (12) 

TTt Usjg A 

When 0 = 7r/2, 



where 


( W l0 2 — ( W 0 + Mg + Oq + ^b) — U Ul (15) 
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and Ql = nA^b^s, u>c/ — are the generally used definitions of the 

Low-hybrid and High-hybrid frequencies (see also (17), (18)). 


In the absence of the eclectic field, i.e when E p = 0, J 2 {p p ) = 1, and 
u> 2 (©) = 0 the three resonance branches are the following: 



— Qg, U>f — U >g, — (Wq + flfl) ~ Wq, 


= 0 , u>l ~ 


Hi /rk . 

, Jl ~ IX — ‘ 1 z jJ 
1 + 7^ 1 + ^f 


W 4 — (^0 + + ^0 + ^ b ) — («I) 2 — 


( 16 ) 


By analyzing the angle dependencies of the dispersion equation D(w,0), 
it is seen that the inequalities Wi(0) < ft B and w 3 (0) < u> B should be fulfilled 
for all the angles 0 < © < 7r/2. However, in [8] the branch u>i(0) given on 
Fig.6 has a maximum at 0 ~ (60 — 70)°, which is comparable with fl p ! What 
is the source of this maximum and of its extremely large value in comparison 
with Cl B is not clear. 


Actually, the numerical calculations of the resonance branch o>i(0) is a 
delicate job because the region of the values of different roots of D(uj, 0) = 0 
is very large. For example, in our case the values of o>i(0) ~ ...10 2 are 
about 10 5 times smaller than u^(0) ~ ...10 7 . Therefore, only by very precise 
numerical calculations this branch can be calculated correctly. As a result, 
we did find that the function u/i(0) is very closely proportional to cos 0 and 
differs from it only in the fourth and following decimal digits. For example, 
at Z = 300 km 
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fig • cos (0.l7r) = 793.466 451 but wi(0.1,7r) = 793.466 349 . 
fig • cos(0.45tt) = 130.513274 but u>i(0.45,tt) = 130.513102. 

and 

D (u>i(0.l7r)) > 0 by wj — 793.46634960 0835 . 

D(u>,(0.1?r)) < 0 by ui = 793.466349600836. 

When E ^ 0, the angle dependence of the branch u>i(0) differs also very much 
both qualitatively and quantitatively from the dependence u>i(0) = f lg — ► 0, 
when E = 0. This is seen from Table I and Figs.l, 2. 

To find the branch i^i(0) close to the cosine rule, when E ^ 0, the 
numerical calculations were done with digits up to the (20-40), and higher 
order. Detailed calculations were also made to track the sharp transition of 
u>i(0) — ► 0, when E p = 0. 

The dependencies of all the resonance branches u>it O 4(0) of the iono- 
sphere for Z = 300 Arm, E ± 0, E p = 0, and p p = 1 and also for comparison, 
the branch a>2(0) for p p = 0.1 are given in Table I and Figs.l and 2. To 
check the revealed here closely proportionality of the branch u>i(0, E 0) 
to cos© shown on Table I adn Figs. 1,2, the calculations were also done by 
the shortened equation (7). Namely, when u 1 -C u >g and cos 2 © m/M, it 
becomes equal to: 

D( u>, 0) = u> 4 u>1(uj 2 l — u 2 ) sin 2 0 — u> 2 ujg(£lg — w 2 )(u;o cos 2 0 — w 2 )+ 
+Q,IuiIu)b(Q 2 b cos 2 0 — u 2 ) [l — -/qC/v)] ‘ cos2 © = 0, 
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Table I. Resonance branches, E =, ± 0, Z = 300Jfcm, 
u 0 = 7.569 • 10 7 Hz, u B = 7.66 • 10®, fi 0 = 7.9 • 10 s , ftp = 5.0857996 • 10 5 , 
Sl B = 834.3, E = 0: w 3 = 7.5694123 • 10 7 , uu = 7.6080677 • 10 7 , p p = 1 


£ 

©° 

E = 0 

u)i,Hz 

E± 0 
u 'i, Hz 

E± 0 


E? 0 
Hz 

E? 0 

u\,Hz 

0 

0 

fts 

ftfl 

ftp 

UB 

>JJ B 

7.5692414 • 10 : 

0.05 

9 

834.299 

824.028 

5.08579 • 10 5 

7.564746 • 10® 

7.564741 • 10® 

7.5701994 • 10 1 

0.1 

18 

- 

793.47 

- 

- 

_ 


0.15 

27 

834.288 

743.37 

5.08575 

6.817977 

6.817917 

7.5772905 

0.2 

36 

- 

674.96 

- 

- 

- 

_ 

0.25 

45 

834.255 

589.94 

5.08563 

5.402849 

5.402727 

7.5886940 

0.3 

54 

- 

490.39 

- 

4.487833 

- 

_ 

0.35 

63 

834.125 

378.76 

5.08516 

3.464215 

3.463922 

7.6000121 

0.4 

72 

833.870 

257.81 

- 

2.357389 

- 


0.45 

81 

832.495 

130.51 

5.07792 

1.119493 

1.194196 

7.6069651 

0.46 

82.8 

831.467 

104.56 

5.07171 

9.58504 • 10 5 

9.57904 • 10 5 

. 

0.48 

86.4 

823.052 

52.38 

4.60609 

4.85080 

5.28379 

7.6077490 

0.49 

88.2 

791.778 

26.20 

2.37225 

2.52242 

. 


0.492 

88.56 

770.529 

20.97 

1.89945 

2.07369 

5.12793 

7.6078754 

0.495 

89.10 

694.912 

13.10 

1.18810 

1.43718 

_ 

_ 

0.499 

89.82 

240.466 

2.62 

2.3772 • 10 4 

8.3068 • 10 4 

5.12217 

7.6078991 

0.4992 

89.856 

195.315 

2.09 

1.9017 

8.1816 

- 

_ 

0.4995 

89.91 

124.165 

1.31 

1.1886 

8.0439 

5.12212 

7.6078994 

0.4999 

89.982 

25.101 

0.26 

2.377 • 10 3 

7.9579 


_ 

0.49992 

89.9856 

20.083 

0.21 

1.901 

7.9566 

5.12210 

7.6078995 

0.49995 

89.991 

12.554 

0.13 

1.188 

7.9522 

- 


0.49999 

89.9982 

2.51 

0.026 

2.377 • 10 2 

7.9544 

5.12210 

7.6078995 

0.5 

90 

0 

10-28 

0 

10-135 

0 

10" 8 

SIl = 

7.9543 • 10 4 

WLp = 

5.12210 • 10 5 

(jJU = 

7.6078995 • 10 7 
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0 ftp pfp 


Resonance branches 
Z = 300 km , P P = 1.0 


E? 0, p p = 0.1 


36° 54* 72° 

Fig 1. Resonance branches of the ionospheric magnetoplasma under the 
action of an alternating electric field $ = E p • cos uj£t , p p ~ E /u>£ 

(see (4)), 0 is the angle between E and the constant Earth’s magnetic 
















and the numerical calculations of Wj(0) become much simpler. The results 
of all these calculations are in a perfect agreement. The value p p = 1 may 
be considered as a strong field characteristic. 


It follows from the figures and the table that the influence of the electric 
field is only large on the dependencies wi(0) and w 2 (©). Quantitatively, the 
behavior of u> 3 (0, E ± 0) differs from w 3 (0, E = 0) only when © > 80° to 
~ 90°, namely w 3 (7r/2, E ± 0) = Sl Lp ~ 5.12 • 10 5 Hz, and u; 3 (tt/2, E = 
0) = fli ~ 7.95 • 10 4 Hz. The influence of the electric field on the high 
frequency branch u; 4 (0) ~ (u>q — ► Uy) is even much smaller. The so called 
hybrid frequency branches u> 3 (0) and u> 4 (0) are also only a few under the 
influence of the ions. Therefore, they can be calculated by the following 
simple equation 


^A — 2 ((^o + w b) ± \J (u>£ +u 2 B ) 2 - 4a^u>| cos 2 0^ , 
which is true for u; 3 (0) only when cos 2 0 m/M. By © = tt/2 
+ n 2 [i — Jq(p p )] + ctg h ^ 


(17) 


u 3 = 


i + S? 


! + 5f 


^ = fiip (18) 


(see (14)). 


What is of a special interest is to find in the experiments the branch 
u> 2 (0), which is of purely electric field origin. Depending on p p , the value 
<^ 2 ( Pp ) is varying from f \ p to 0. By small values of E p or large frequencies 
u>e, when p p << 1, and «/ 2 (p p ) ~ (1 — p 2 / 2), even ELF parametric waves 
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o> 2 (0) may be excited. Therefore, by theoretical calculations of u> 2 (0), the 
thermal velocity of the particles should be taken into account, and the kinetic 
approximation should be used. It brings to the following equations: 


U) 2 = COS 


e( noft-JoM + 3tyV (19) 

V 1 + » A ■ sm! © ' [1 - JoMl / 


JSM1 

and the Landau damping of these waves is equal to 


7 = 


x ur 


8 I*# 


h' Ezp 


2 kVj 


+ cos 2 0 - 


W 4 u3 J ’ 


(20) 


when kxPd < 1 (see the Table on p.4), v, <C u*|| < v e , taking also into 
account that D e k e <C 1. By p p > 1, even when p p ~ 1, 

n 2 cos 2 0 


u 2 ~ 


+ 3k 2 v 2 cos 2 0. 


1 + sin 2 0 
When -C 1 and 0 = 0 (isotropic plasma), 


( 21 ) 




= n? = n2 (i - •£(»,)) +(*»0 2 = i • ^’•(t-4) 2 +3(to i ) 2 - ( 22 ) 


This equation, rewritten as 

n; = 3tx + *’?), (23) 

demonstrates a very important physical result. Longitudinal ELF waves, 
propagating along the direction of the electric field E p , can be parametri- 
cally excited in a plasma. Their velocity V p is determined by (23). To find 
experimentally this branch of ELF parametric waves in the ionosphere, it 
would also be a great achievement of any satellite mission. 
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Results of calculations of the resonance branches u; 2 (0) and cj 3 (0) at 
different altitudes Z of the ionosphere and different values of p p ~ E p juj\ are 
given in the Tables II and Ila. The characteristic frequencies u> 2 (0 = 0) = fi p 
and w 3 (0 = 7 t/ 2) = ft£, p of these branches are of pure electric field origin, 
they are given in the following table: 


Z, km 

PP ~ 

0.1 

Pp " 

0.5 


ft p 

f^Lp 

ftp 

ft^p 

200 

1.89 • 10 4 

5.60 • 10 4 

9.25 • 10 4 

1.05 • 10 5 

300 

5.58 • 10 4 

9.70 • 10 4 

2.73 • 10 5 

2.83 • 10 s 

400 

CT> 

b* 

00 

t— ‘ 

o 

•fr. 

1.18 ■ 10 5 

3.21 • 10 5 

3.34 • 10 5 


Z, km 

Pp - 

1.0 

Pp ~ 

1.5 


ftp 

ft/,p 

ftp 

ftip 

200 

1.72 • 10 5 

1.77 • 10 5 


E 

300 

5.08 • 10 5 

5.12- 10 s 



400 

5.99 • 10 5 

6.04 ■ 10 5 




To illustrate the dependence of w 2 (© = 0) on p p (see (22)), the ratio of 
ftp/fto is also given in the following Table lib and on Fig.2a. 
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Table II. Resonance branches, E ^ 0, 0° = 


n in 
^ - 

U)b = 
8.02 • 10® 
7.66 • 10* 
7.31 

7.61225 

7.28151 

6.94849 

6.44321 

6.18611 

5.90250 

4.65212 

4.48749 

4.28120 

2.43419 

2.35652 

2.24813 

1.95781 

1.89634 

1.80926 

1.47453 

1.42902 

1.36377 

N o 

^ r-J 

% % 

II ~ 

m ’ • « 

3 § 8 «-• 
oo r-' 

7.61225 

7.28152 

6.94851 

6.44234 

6.18617 

5.90260 

4.65218 

4.48764 

4.28144 

2.43436 

2.35692 

2.24877 

1.95803 

1.89686 

1.81008 

1.47483 

1.42973 

1.36488 

Hz 

0.5 

CO CO 
M ° ° 

~ • • CO 

3 © <o 

00 N 

7.61226 

7.28154 

6.94853 

6.44236 

6.18622 

5.90268 

4.65224 

4.48778 

4.28165 

2.43452 

2.35572 

2.24929 

1.95822 

1.89728 

1.81071 

1.47509 

1.43025 

1.36558 

N 

SC! o 
- II 

3 d. 


7.61226 

7.28154 

6.94854 

6.44237 

6.18624 

5.90271 

4.65226 

4.48783 

4.28173 

2.43457 

2.35738 

2.24980 

1.95830 

1.89743 

1.81093 

1.47519 

1.43042 

1.36580 

Hz 

1.5 

ftp — 

2.30238 • 10 5 
6.78683 
7.99807 

2.30239 

6.78681 

7.99805 

2.30238 

6.78677 

7.99798 

2.30236 

6.78664 

7.99770 

2.30226 

6.78583 

7.99602 

2.30218 

6.78515 

7.99450 

2.30200 

6.78340 

7.99012 

0T 

Z H 

to 

O 

H O fli 

H ^ 8 8 

AM °0 Oi 
C O Oi 

«5 w 

?“H 

1.72534 

5.08578 

5.99336 

1.72533 

5.08571 

5.99323 

1.72530 

5.08549 

5.99278 

1.72515 

5.08419 

5.99016 

1.72499 

5.08316 

5.98798 

1.72471 

5.08071 

5.98258 

Hz 

0.5 

ftp s 

9.2532 • 10 4 
2.7276 • 10 5 
3.21429 

9.2533 

2.72757 

3.21427 

9.2532 

2.72751 

3.21416 

9.2530 

2.72733 

3.21382 

9.25160 

2.72633 

3.21183 

9.2505 

2.72566 

3.21030 

9.24830 

2.72386 

3.20686 

tq ® 
. II 
3 a 

ft P = 

1.8938 • 10 4 
5.5823 
6.5784 

1.8938 

5.5824 

6.5785 

1.8939 

5.5824 

6.5786 

1.8940 

5.5822 

6.5782 

© © 
w o ^ 

05 OO K 
00 tO iC 

H IO <0 

to fO N 
CO 00 O 
C 7> t>- 
00 to to 

H to © 

O (O Tf 

co ^ eo 

05 N ^ 

00 to to 

i—i to CO 

1 

N 

200 

300 

400 

3 O © < 
3 © © < 
N M ^ ( 

§88! 
N rt ^ ( 

3 0 0 
3 O © 

N rt ^ l 

3 O © < 
3 O © 

N « ^ ( 

§ § 8 
N rt ^ 

8 8 8 
N M ^ 


3 

oro 

0.20 

0.30 

0.40 

0.42 

0.44 


17 









































Table Ha. Resonance branches, E ^ 0, 0" = ir/3 (Continuation) 


Hz 

1.5 

9.86129 • 
9.56961 
9.1684: 

4.9467! 

6.8108 

8.02121 

00 N ^ 
00 0*0 
TT O ^ 
115 00 O 

CM CO 00 

2.3183! 

6.79951 

6.04011 

2.3182! 

6.7990- 

2.3182 

6.7990: 

6.040 H 

II CM © — 
ft. 00 2 ° 
^ CO N o 

^ cm cd cd 

N O 

st: ~ 

9.86584 • 10 5 
9.57904 
9.16468 

4.95610 

5.28379 

6.09864 

2.53076 

5.13210 

6.04777 

1.77297 

5.12217 

6.04016 

1.77228 

5.12215 

6.04011 

1.77207 

5.12210 

6.04010 

flip = 
1.77206 
5.12210 
6.04009 

so 

Z H 

9.86981 • 10 8 
9.58369 
9.16371 

4.96299 

4.86919 

4.73267 

2.52652 

2.99643 

3.44098 

1.05652 

2.82860 

3.34378 

1.05140 

2.82807 

3.34334 

1.04983 

2.82790 

3.34320 

ft Lp = 
1.04977 
2.82789 
2.34319 

st : ° 

. ii 

S £ 

9.87130 • 10 8 
9.58499 
9.16349 

4.96534 

4.85135 

4.66967 

2.52561 

2.52816 

2.49773 

— 0 

2 2 

o 

(v ^ 

© 04 Oi 

. T— 1 

2 2 

o . 

N J2 J 

CT5 t~- 

£ oo 

« 2 

’2 S 

’ n nr, 

H O S 
f- O 

g ® s 
• 1— * 
•» -4 

_ >e 

2 2 
II 7 2 • 

SS 2 g 

id ^ 

** U5 

Etl -J 

O 

H ic N 
* © CO 

00 2! 22 
© 

£ » ^ 
cm 

2.29801 

4.76855 

4.55145 

2.23096 

2.38895 

2.27853 

2.4351 • 10 4 
2.3899 
2.266 

1.2266 

1.195 

1.1330 

*"• © © 
• © © 
eo co cm 

^ cs c4 

CM 

© 

Hz 

1.0 

ie 

O 

H H w 

• r- oo 

N 00 

oo to 

n o a; 
id id 

1—1 

1.72883 

4.60609 

4.48584 

1.68376 

2.37225 

2.26287 

** 

o 

^ £5 

£ « « 
o ci ^ 

CM 

1.2023 

1.1886 

1.1330 

^ t- © 
• © 
to CO CM 
o <N CS 

CM 

© 

SO 

Z H 

A 

2 o -i 

. t- 1-H 

OO CO 

22 ^ 

^ 

(N <N W 
© 

9.2056 

2.6806 

3.10017 

9.0454 

2.17906 

2.13298 

2.1634 ■ 10 4 
2.3087 
2.1953 

05 tO 00 
CO T** t* 
00 to 05 

© ^ © 
1-H »— H 

n 

© 

^ © © 
• © © 
r— co »— • 

q « 

CM 

© 

* 

&3 ° 

~ II 

a. 

d ^ 

r 

O 

^ O t*- 
• Tf b- 
iO © ^ 
h H5 W 

§o « <© 

rH 

1.8835 

5.5070 

6.4317 

« « ^ 
cm co *o 
to 00 ^ 
00 CM © 

id cd 

o ' 

2 ^ CO 

• OJ 
§ O N 

& 

H ^ N 

* to to 

o 00 CO 

§ <0 cd 

"O 4 

o ^ 

5 ^ *o 
• * I s " 

^ ^ CM 

© K. • 

« ^ 

© 

i ! 

HI 

■ 

■ 




200 

300 

400 


1.46 

00 

rr 

© 

0.49 

0.499 

0.4995 

0.4999 

SO 


18 















































Table lib. Frequencies of the longitudinal component (0 = 0) 
of the parametric resonance branch w 2 (0) 


ftp _ w 2 (0 = 0) 
fto fto 


r 72 / xI 1 / 2 27r e E p -= _ 

[1 - J 0 (pp) j , Pp = , E = £/ p COSWE(f), WE > Wo 


Pp 

0 

0.02 

im 

0.1 

0.5 

1.0 

mm 

2.0 

MPp) 

l 

0.9999 


0.9975 



0.5118 

0.2138 

£lp/£l 0 

0 

0.01 

0.024 

0.071 

0.345 



0.976 

Pp 

2.4048 

3 

3.832 

5.0 

5.520 

6,0 

6.5 

7.016 

MPp) 

1.0 

-0.2601 

-0.4028 

QSS1 


0.1506 

0.2601 

0.3001 

n p /n 0 

max 

1.0 

0.966 

min 

0.915 


max 

1.0 

0.988 

0.965 

min 

0.954 

mm 



8.5 

8.650 



10.0 

10.713 

Jo{pp) 

0.2663 


0.0419 

0 

-0.0903 

-0.1939 

-0.2459 

- 0.2497 

Slp/Sl o 

0.964 

0.985 

0.999 

max 

1.0 

0.996 

0.981 

0.970 

min 

0.966 


III- 2. Parametric VLF resonance instabilities 
and oscillations. 


The principal properties of the parametric resonance instabilities are, 
more or less, similar in an isotropic ( B = 0), and anisotropic ( B ^ 0) plasma. 
In an isotropic plasma the resonance regions of the parametric effects are of 
a special interest when we ~ sw 0 , s = 1,2,.., i.e. when the frequency of 
the electric field is comparable with the electron Langmuir frequency or its 
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Fig 2a. Dependence of fl p /fio of the longitudinal component 
of 0^(0 = 0 ) = fl p on p p . 
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harmonics. In an anisotropic plasma, distinguish resonance instabilities arise 
when the overtones of the electric field 

sue — w 3 (0), ~u> 4 (0), ~ [w 3 (0) + w 4 (0)], (24) 

i.e. by approaching of the resonance hybrid frequencies (see (17), (18)) or 
of their sum. However, the calculations of the resonances in a magneto- 
plasma are considerable more complicated. And for the ionosphere detailed 
numerical calculations are given here only when w# ~ su o, i.e. not taking 
into account the influence of the magnetic field. These results give a general 
qualitative and, in a sense, also a quantitative presentation of the behavior of 
the parametric resonance instabilities. Some equations of these instabilities 
in a magnetoplasma are discussed for illustration at the end of this section. 


a. Isotropic plasma ( B = 0). The high frequency spectrum of the dis- 
cussed here parametric oscillations has no peculiarities. When the frequency 
of the electric field is not much larger than the electron Langmuir frequen- 
cy u o, we can approximately use for the calculation of the spectrum of this 
branch of oscillations the equation 


u 




)■ 


(25) 


It becomes 


when u>g >• Wq. 


U — + Jq(Pp) ' 


(26) 
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However, the low frequency band w < w 0 of these oscillations becomes of 
resonance nature 


-oo ( w + SUJ E) 2 ~ 


(27) 


(see [8]). Let us call the region, which is far enough from the resonance 
region, the out of resonance region - Out. The equation (27) for this Out 
region becomes 


'f Jim ■ , , = *5 (i - *,(*» , 


(28) 


7 TX 


where x = u>o/u>e, and 

$x{Pp) ■ Jx (Pp) ' J— xi^Pp)- (29) 

sin x 

In (25) to (29) are Bessel functions. Equation (28) was used here for 

numerical calculations of u>/fio and 7/fio of the oscillations u 7, and growth 
rates 7 of these waves when u> 2 becomes negative. By the converse case, 
namely when 

((w + sue) 2 - Ho) < w l , 

let us call it the resonance region - Res, the solution of (27) brings to the two 
following formulas which estimate the oscillations u and the growth rates 7. 
Namely for s = 1, when ujq/ue < 1, 

Meff 


— SL _ l J Mill 

Ho Hq 4 


»/*l 1/3 


m 


|±(x J - l) 3 + (32 . ■ |(x 3 - 1)|)‘ /2 | 


(30) 


and for ojq/uje > 1, 


+ (<** - »* + 32 — ■ *<*> ■ k* 2 - di) 


1/2 


( 31 ) 
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where M e fj is the effective mass of the ions of the plasma and the signs 
< ± > in the {...} brackets were used by calculation of w/flo and 7 /^ 0 - 
Certainly, the frequency region J- 1451 - — 1 1 of applicability of the Out and Res 
formulas are not known precisely. Approximately, they might be estimated 
by the crossing of the Out and Res branches which are seen on Fig.5 (see 
below). 

Results of numerical calculations by equations (28) of the dependencies 
w/flo and 7 /wo in the Out region, and by equations (30), (31) in the Res 
region for different values of p p are given on Figs. 3 and 4. The values of the 
growth rates in the resonance region are specially selected in Table III, and 
in details they are given for the Out resonance region in Table IV . It is seen 
that for all of the values of p v used here, the maxima of the growth rates 7 
axe comparable with the ion Langmuir frequency. Both dependencies are also 
presented for the second harmonic, s = 2 on Fig.5 for p p = 1. The intersection 
regions of the Res and of the Out dependencies give a rough presentation of 
the values (wo/w#) of applicability of the Out and Res formulas (28), (29) 
and (30, 31). They are equal to 

~ (0.9 -1.1), (— ) ~ (1.1 -1.8). (32) 

\UE/R e , WE /Out 

It was very important to evaluate even roughly the amplitudes E p of the 
electric field which may generate the resonance oscillations. For this purpose 
we used the asymptotic solution of the dispersion equation in the kinetic 
approximation for an isotropic plasma. If u>o / u>e < 1, the threshold value of 
the amplitude of the electric field E Pttfire ,h, which is necessary for excitation 
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Fig 3. Dependencies of the oscillations u> and of the growth rates 7 , 
Pp ~ EpI^E, u 0 is the Langmuir frequency, u>e is the frequency 
of the electric field. 
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Resonance region- 
- Res. 



Pv = 1.5 



Table HI. Resonance regions - (Res.) 


" = ( “ ,/4) '(^) '(- (l3 " 1)2 + ■ jj 

*<»=!: aj = 1, 02 = 0; x >, = 1 : Oj = y/2, o 2 — 1 


908 

920 

935 


0.945 

0.954 


maxi :0.9632 


0.9690 


0.977 


0.9794 


0.982 

0.984 

0.986 

0.991 


0.9929 


1.00 

1.002 

1.005 


1.0088 


1.009 

1.015 


1.026 


1.039 


max2: 1.046 


1.10 

1.15 

1.20 


p p = 1.5 


03 

80 

09 


1.20 

1.26 


1.291 


1.882 


1.161 


0.04 

0.80 

0.93 

1.02 


1.08 


maxl:1.105 


1.09 


max2: 1.607 


1.56 


25 

01 


0.94 


maxl:0.7423 


0.851 


.980 

.046 


max2: 1.080 


1.065 


1.01 


0.70 

0.60 

0.52 


009 

.14 

.18 

.25 


maxl:0.262 


0 

.30 

.36 


max2:0.376 


0.37 

0.35 


0.29 


0.24 


0.22 


[i] 



























































Table IV. Resonance outside regions - (Out.) 


y = £■ = (i - gjfr • Up,) ■ U(p,)) 1,3 , p, = — ■£$■■ y = 2t y 

Uq \ sinirx / m wi A A 


Pp — 0.1 


w/fl 




jy/%. 

0.247 

0.244 

0.235 

0.221 

0.202 

0.196 

7/fto 

0.192 

0.160 


0.656 

0.673 

0.699 

0.738 

0.797 

0.889 

1.052 

1.195 


0.874 

0.896 

0.928 

0.976 

1.048 

1.162 

1.366 

1.544 


1.440 1.850 

2.002 2.556 

3.133 3.984 

4.416 5.607 

7/fto 
48826 
— ► oo 
5.551 
3.906 
2.432 
1.673 
1.112 
0.840 
0.491 
0.405 
0.298 
0.072 

ui/fto 


.052 
|.138 
.264 
.412 
0.181 I 0.549 




0.038 0.639 

O>/ft 0 u/fto 

0.667 

0.104 0.673 

0.352 


0.062 0.631 



1.96 

1.97 
199 
1.999 

1.999 

1.9999 

2.00 


2.01 

2.10 

2.70 

2.90 

2.93 

2.935 


2.938 

2.94 

2.96 

2.98 
0.2983 

2.985 

2.99 
2.995 

2.999 

2.9999 

2.99999 

3.000 


p p = 0.1 0.5 


7/Qq a; /A 

9 0.077 

3 0.149 

0.385 


1.0 


U>/A 


Jt/Ho 

0.0250 

0.0249 

0.021 

7/flo 


1.52 
— ► oo 
0.025 
0.0245 
0.0241 
0.0239 
0.0236 
0.0232 



w/flo w/A 0 
1.354 5.127 

4.324 

26.78 10055 

— ► oo 

7/flo 
1.661 
0.608 
0.299 
0.253 
0.240 


0.134 

0.126 

0.123 

0.0122 

0.121 

0.118 


w/A 


1.586 

1.848 

3.256 

w/Ao 

10.371 

20310 
—*■ oo 
7/A0 
3.307 
1.104 
0.436 
0.260 
0.120 
0.034 


7/A0 I u//A 



0.428 

0.298 

0.274 

0.268 

0.262 

0.253 


u*/Ap 
5335 
— *• oo 
1.129 
0.614 
0.508 
0.487 
0.463 
0.440 
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Fig 5. Dependencies of the frequencies u> and growth rates 7 in the 
resonance - Res and out-resonance - Out regions for p p — 1. 
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of this branch of oscillations, is determined by the formula 


= — (4tJ NkT, + 4xJV*7i) = Hfii . NkT',, (33) 

U^o 

and the threshold amplitude of the electron oscillations (see (4)) is deter- 
mined by 

r , = ± S} 

E^thresh Q^AT^'V 


87 tNkTc 

In the ionosphere the values of E p<t hresh and r E<thTeah are equal to: 


Z, km 

B? hr .. h CGSE , 

E V/m 

r hhr e ,h r E cm 

200 

3.47- 10" 10 , 

0.56 

1.51 • 10“ 5 , 3.9 

300 

05 

1 

O 

iH 

00 

Tt; 

rH 

1.15 

1.09 10~ 5 , 3.3 

400 

8.13- IQ' 10 , 

0.86 

1.42 -10" 5 , 3.8 


It is seen that E Pyt hresh ~ 1 V/m. The used by these calculation values of 
* 'eii 4ttA^ 7;, and D e are given above in the Table of Section II. The region 
of the applicability of (33) an (34), 

0 < u E — \Ju>l + n 2 0 < \k 2 D\ (35) 

is very narrow and depends on the wave number k = 2ir/X. 

When u>e < \Ju > 0 + Hq the threshold of the electric field is equal to 
E p , thresh = 2 • (A ■ 4wN.kT, + 4tt N t KTi) 1/2 = 2\/2 • (A • 4irN e KT e ) 1/2 , (36) 


29 




where A = ^1 — uje !\fw o + ^oj- I- e - 


J p } thresh 


N 


\Jwq + — we 

\j' wi + 



( 37 ) 


and it is increases with ujo/we- In the regions of the maximum value of the 
growth rates 7/flo of the VLF resonance oscillations (see Figs. 3, 5, and 
Table III) the values of E Vtt hresh are equal to: 


Z, km 

CGSE 

Ep } thre»h V/m 

200 

5.46 • 10 4 

16.38 

300 

1.27 -10" 3 

38.1 

400 

1.27- 10" 3 

38.1 


They are much larger than the values estimated by (33), and smaller than 
the critical values of E pfi given in the table of Section II (see (3)). 

Thus, the amplitude of the electrical field which can excite the parametric 
resonance branches u>\ t0 4 i analyzed in Section III-l, and the VLF resonance 
oscillations, considered in this section, are changing in different conditions 
from a fraction of 1 Vjm to (16 — 40) V/m, and may be equal to some 10’s 
and even more Volts per meter. 

The maximum growth rates of these oscillations described by (33) and 
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(34) are the following; 


wo < l /_7_\ _ M ' (rh ~ r^thrcsh) ~ M r\ 

VE ’ \ / max rn r 2 De + r 2 Di “m 2^’ 

where the amplitude of the oscillations of the electrons is: 


(38) 


® Epjhresh 
T E = 2 

m u>£ 


When u>o/u>e > 1, than 

( 2 _\ g> »gy /3 

V ^0 / max ^0 \ 2 U}q ) 


In the ionosphere the values of {l/^o) max are the following: 


(39) 


(40) 


Z, km 

cm 2 

r E,thre*h 1 c ™ 2 

( ll^o)max 

Mq/^E < 1 

UJq/u>e > 1 

Pp — 1.5 

© 

II 

p p = 0.5 

Pp = 0.1 

200 

6.14 • 10 ~ 3 

1.51 • 10 " 5 

0.34 

2.85 

2.43 

1.64 

0.56 

300 

3.514 • 10" 3 

1.05 • 10 ~ 5 

0.61 

2.56 

2.18 

1.46 

0.50 

400 

3.151 • 10 ~ 3 

1.42- 10 " 5 

0.20 

2.27 

1.93 

1.29 

0.44 


By calculating ( 7 /fto) mox we used the table given in Section II. It is seen 
that the growth rates 7 max are rather large, and are comparable with the 
ion Langmuir frequency fl 0 - I.e. the VLF parametric resonances in the iono- 
sphere are phenomena of great vitality. 


Let us note here that the formula (40) of ( 7 /fIo) max is of pure cold plasma 
approximation origin. The values of the growth rates, given in the table are a 
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little smaller than on Fig. 3 because they were calculated by using the mass of 
protons (see Section II). Vice-versa, when ujo/SIe < 1, they are smaller than 
the values given in Fig.3. However, the data given here are namely the result- 
s of the more correct kinetic approximation, which brings to the formula (38). 


b. Anisotropic plasma (B ± 0). When the frequency of the electric 
field su>e , s = 1, 2, .. is not much larger than the hybrid resonance frequencies 
u> 3 (0), u; 4 (0), and cos 0 >• m/M, the spectrum of this low frequency band 
w<«£ may be described by the equation 


<4 = 


nj, + ng-*(a; g ,fc)] 


yj • 4>(w b , £)] 2 - 4f^fig$(u; iJ , £) cos 2 0 


(41) 


where 


$(u E ,k)= ^2 J?(p p b) 


s 2 u 2 e (s 2 u: 2 e - u%) 


(42) 


(sW E - - u, 2 ) 

and a> 3 (0) and w 4 (0) are estimated by (17). Equation (42) is identically 
equal to 

$(u>E,k) = 1 -| 2 [^*ii (PpB ) — $* 2 (Ppb)} ~ ^ x 2 (Ppb)i (43) 

u> 4 

(see [8]), where 


,i{Pvb) — 


sin xxi t 2 


( PpB ) {PpB)i 


and xi = a ; 3 (0)/w£, x 2 = u 4 (Q)/ue. 


(44) 
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Under the influence of the magnetic field ~B the characteristic parameter 
of the plasma (see (4)) becomes equal to 

= (45) 

where 

sin x sin © sin <p) 2 
(46) 

and 0 is the angle between the magnetic field B and wave vector fc, x is 

— * — t 

the angle between B and the electric field E, and tp is the angle between the 
surfaces ( k,B ) and (E P ,B). 


*(x,V?>9) = 


cos X cos 0 + 


u>£ sin x sin © cos tp 

( W E — ) 




I 


From (41) and (42) it is seen that when sw 3 < <*>3(0) and 0/4(0 ), $(u>e, k ) 
and the roots of (41) become negative, namely u> 2 (0) and a^(0) < 0. I.e. 
instabilities arise in the magnetoplasma and, particularly, the interested us 
here oscillations of the hybrid band 013(0) grow up. Their growth rates are 
maximal when ps — 1 and become close to the roots of hybrid branch 0/3 
(see (41)). Namely, when o/ 2 


2 _ SU} E l 4 f 2 o(s 2 o>g — ^ b\J^{Ppb) 


= 


s 2 u: 2 e {2s 2 u}% - <x>l - u>b 


2/3 


a±, 


(47) 


where a± is a numerical coefficient ~ (1 to 2) (see [8]). 


It is important to note here that because of the dependence of the hy- 
brid resonance frequency on the angle 0, the width of the resonance region 
in a magnetoplasma is very much broader than in the width of the reso- 
nance regions in the ionosphere (isotropic plasma) calculated above. This 
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circumstance facilitates the discovery of these resonances in the ionosphere 
by experiments in Situ on satellites. 


In the threshold region of the instabilities in the low frequency band 
\ u + *7 1 u*(0), when Ao> 3 = (u> E - u> 3 ) < 0, the threshold value of 

the amplitude of the electric field E Pit hreshi which is necessary for excitation 
of this branch of oscillations, is equal to 


F 2 __ 

^p, thresh 


- • NkT» • 


Wn 


1~4 

1 + 2 4 

i-4 


and the maximal growth rate of these oscillations 

w E u>% sin 2 © 


v- 


7mo * - y I eI 


Zl 


- 1 


1 + 


~'p,threah / \ (w| - ) COS 2 0 + U/ 

(see the formulae (33) to (38) given above). 


E sin 2 0 ) ’ 


(48) 


(49) 


Summary 

In this report results of numerical theoretical calculations of e.m. os- 
cillations and waves produced in the ionosphere under the influence of the 
electric field E — E p cos uj£ t are presented. They are chiefly done for the 
VLF - (Hb < w — + fi 0 ) and ELF - (0 < u> < VIb) frequency bands, where the 
parametric activity of the electric field is especially strong (fig and fi 0 are 
the gyro- and Langmuir frequencies of the ions). 
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The resonance branch of pure electric field origin u>2(0) = (Q p to 0 is of 
a special interest. The frequency fl p = fi 0 (l - Jl{p v )) is changing, depend- 
ing on the amplitude E p and frequency u>e of the electric field, from to 
zero. I.e. parametric longitudinal ELF Alfven waves can be generated in the 
ionosphere. Another important result is the found dependence close to the 
cosine of the ELF resonance branch u>i(0) = (f Is — » 0). This is in a strong 
contrast to the angle dependence of u>j(0), when the electric field is absent. 
Therefore, this effect may be an indicator of the existence of the electric field 
in the ionosphere. The necessary amplitudes E p of the electric field to pro- 
duce the effects discussed here are changing in different cases between about 
1 V/m and some 10’s of Volts per meter, and more. 

To discover and to investigate these effects by satellite experiments will 
be an important contribution to the ionosphere knowledge and to plasma 
physics. 

Acknowledgement : I am very grateful to Svetlana Alpert for the excep- 
tionally helpful and laborious computer calculations. 
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I . INTRODUCTION 


The ionospheric measurements of the BICEPS observation plan have a scientific 
relevance in themselves, in terms of improved knowledge that through them will 
be gained on the ionosphere and the magnetosphere. 


Additionally, these measurements have an important application potential to communi- 
cations, could the BICEPS measurements be correlated ( and this is a distinct possi- 
bility) with the degradation of the communication waveform (due to ionospheric effects) 
experienced by satellite microwave links. These microwave paths extend from geo- 
synchronous height to ground and do cross ( at the same time that the BICEPS measure- 
ments are made) the volume of the ionosphere where BICEPS operates. 

With the ever-increasing bandwidth requirements in satellite-to-ground microwave 
communications, and with the continuous progress in communications equipment techno- 
logy (that makes gigahertz-bandwiths feasible) , the frequency-spread and time-spread 
characteristics of transionospheric propagation paths are fast becoming the limiting 
factors in achievable communications performance, in terms of data rate and error rate. 

An ionospheric science mission such as BICEPS, if coordinated with existing broad-band 
microwave links from synchronous heights to ground, would make a vital contribution to 
the understanding of the limitations imposed by the ionosphere on microwave link per- 
formance. Ionospheric effects would be measured in real time , for comparison with their 
influence on the waveforms transmitted at microwave from synchronous heights to ground 
(with their propagation paths intersecting the ionospheric regions where BICEPS performs 
its measurements) . 


Once that causes and effects are well understood , it will be possible to design 
a signal processor, to be added to the ground terminals of the satellite microwave 
H n k, with the goal of reducing as much as feasible the limitations imposed by the 
ionosphere on link's performance. 

The channel is time spread and frequency spread , and a full correction would require 
a "phase conjugation" approach, never tried at microwave . However , it might be possible 
to correct for gross features, such as the waveform degradation experienced by a short 
pulse. An inverse filtering approach would be quite effective in correcting for this 
degradation, we would expect. 


II ~ BICEPS Ionospheric Measurements 


The basic experimental arrangement for INVESTIGATION B is a multifrequency phase- 
coherent link established between the two terminal of the pair of BICEPS satellites. 
^Id^S ( Bistatic Canadian Experiment on Plasmas in Space) could use two separate free— 
flyers called WISPRS ( Wave in Space Plasma Research Satellite) or the tethered pair 
BOLAS (Bistatic Observations using Low-Altitude Satellites). For WISPRS, separation 
distance is from a few meters to several tens kilometers, while for BOLAS is, at most, 

1 kilometer (the length of the tether). However, the tether will be cut at the end of* 
the initial series of measurements, and the two satellites will drift apart freely up 
to large distances, making the configuration as suitable to general-purpose bistatic 
measurements as WISPRS. 
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The scientific objective of INVESTIGATION B is to detect and characterize electron 
density perturbations in the ionosphere, such as large-scale to small-scale dis- 
turbances, T.I. D.s, irregularities producing scintillation, spread-F phenomena, 
etc. The frequency band of the bistatic sounding will be 0.1 MHz to 20 MHz. 

Specific phenomena to be obserbed are the following. 


Travelling Ionospheric Disturbances an d Atmospheric Gravity Waves 
gravity waves (AGWs) are an important mechanism by which energy and momentum are 
transported from one region of the atmosphere to another. Their ionospheric manifestations 
are called Travelling Ionospheric Disturbances (TIDs) [e.g Hines, 1974]. There are many 
sources of AGWs: auroral, meteorological, stratospheric and mesospheric heating, and 
geological. In the thermosphere the energy carried by these waves can be comparable to the 
input from solar EUV [e.g. Ching and Chiu, 1973]. AGWs may grow non-linearly, 
generating smaller scale structures in ionization and turbulence in the atmosphere re. 2 . Cos ta 
and Kelley, 1978]. 


One of the most important sources of thermospheric TIDs is Joule heating by auroral 
precipitation, generating wavelengths from tens of centimetres (acoustic waves) to thousands 
of kilometres (gravity waves). These waves propagate to mid-latitudes and are believed to be 
a cause of mid-latitude spread F [e.g. Bowman, 1981]. The direction of arrival 
measurements provided by the Doppler capability [Dyson, 1978] of the BICEPS receiver will 
provide a new ability to map details of Til) structure from satellites and relate them to 
sources and spread F. AGWs dissipate as they propagate to hi gh altitude and their 
interaction with the ionosphere causes the effects of TIDs to be greatest near the F-layer 
peak. Tins generally have horizontal wavelengths greater than 100 km but structure with 
smaller wavelengths occurs in association with spread F. Consequently the BICEPS TED 
study will be capable of observing horizontal structures spanning tens to thousands of 
kilometres. 


The specific aims of BICEPS TED study will be to: 

• Determine TTD wave characteristics near the F-layer peak (wavenumber, direction 
of travel, frequency and amplitude); 

• Locate sources of the waves; 

• Infer the energy and momentum transferred by TIDs and the AGWs driving them; 

• Determine the occurrence and structure of s mailer-scale disturbances genera ted bv 
TIDs. 

An important aspect of achieving these aims will be simultaneous observations by 
ground-based ionospheric radars. 
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Mid-latitudf F Generally raid-latitude spread F is not as severe as that 
occurring at low and high latitudes and there is still no accepted theory explaining its cause. 
Many measurements indicate that the structures responsible are "ripples" with horizontal 
scale-sizes of tens of kilometres [e.g. Bowman, 1960; Scali, 1989]. However other 
measurements indicate smaller irregularities may sometimes be significant [Bowman etal, ^ 
1988; Dyson and Roberts, 1989], The situation is further complicated by occurrences of very 
severe spread F but these are clearly associated with irregularities moving equatorward from 
higher latitudes [Scali and Dyson, 1992]. BICEPS will offer a unique opportunity to study 
mid-latitude spread F. The BICEPS sounder will detect spread F and the measurement of the 
direction of arrival using the Doppler shift will allow us to distinguish "ripple" structure. The 
latitude coverage will enable the boundary of high latitude spread F to be determined, as 
will the analysis of the amplitude statistics [Scali and Dyson, 1992]. The Langmuir probes 
/-in be used to dete rmine the irregularity spectrum, crucial information in determining 
whether plasma instabilities are involved. Observations near perigee are prefer red, and 
co-ordinated ground-based ionospheric radar measurements will be important 


t arge-scflle structure Complicated large-scale ionospheric structures occur in the 
transition from typical mid- to typical high latitude conditions. At night the transition is 
marked by a pronounced trough in electron density and minor troughs can occur as well 
[Muldrew, 1965]. Total-electron-content (TEC) [Davies, 1969] observations reveal large- 
scale ionospheric structure and recent measurements over southern Australia show an 
fowj ra firm enhancement in the early evening possibly related to equatorward movement or 
expansion of high latitude convective circulation [Breed et al., 1992]. In early winter spread 
F at mid-latitude stations can be primarily due to equatorward movement of high latitude 
ionospheric structures [Scali and Dyson, 1992]. 

The BICEPS sounder will be used to study the features of the mid- to high-latitude 
transition that are not understood. The sounder has the advantage over other measurement 
techniques of being able to measure both horizontal and height variations in the ionosphere. 
The Doppler-direction of arrival capability will enable vertical echoes to be distinguished in 
the presence of spread F, further enhancing the ability to determine large-scale ionospheric 
structure. Simultaneous measurements with the^ Langmuir probes will provide measurements 
of electron temperature, enabling regions of ionospheric heating to be identified. The particle 
detectors will also identify regions of energy input. 

Coiar control phenomena A particular experiment that capitalizes on the bistatic 
configuration is the measurement of the refractive Doppler shifts (using phase-coherent 
carriers of the HF link) caused not only by well-known large scale features of the ionosphere 
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(such as the equatorial geomagnetic anomaly), but also by lesser-known features, such as the 
sharp peak in electron density that occasionally appears at the day-to-night transition, and 
that may be due to terminator-related shock waves. The latter features were observed in 
some of the orbits of the Apollo-Soyuz Test Project (ASTP) spacecraft-to-spacecraft 
ionospheric sounding [Estes and Grossi, 1984], 

Figure i shows an example of the results of the ASTP space-to-space sounding. 
Hie day/night transition and magnetic equator crossings are indicated in the figures. This 
marking is almost superfluous for the day/night transitions, since the change in election 
density is so dramatic. Daytime crossings of the equator are from south to north, and 
nighttime crossings from north to south. 

A prominent feature seen repeatedly in different revolutions is the geomagnetic 
anomaly or equatorial trough first observed by Appleton [1946] in which the electron- 
depleted magnetic equator is flanked by crests to the north and south. In addition to the well- 
known features mentioned above, i.e., day/night ionization differences and the Appleton 
anomaly, there is another feature that should be noted. In the two consecutive revolutions 
displayed in Figure l there is a sharp peak in electron density just on the day side of the 
day-to-night transition point. The peak is especially prominent in the second revolution. 

This phenomenon may be related to the motion of the terminator with respect to the 
ionosphere. The proposed high-inclination orbit and telemetry facility for BICEPS mean that 
the solar terminator will be routinely observable once per orbit on a significant number of 
occasions throughout the mission. 

It has been suggested that a shock wave generated by the supersonic motion of the 
terminator could be a source of HD, or of other irregularities associated with spiead-F 
phenomena. The condition for supersonic motion to occur is [Cole, 1974]: 

> K 

where V T is the speed of the terminator with respect to the rotating earth coordinate frame, 
is the wind speed (positive toward the east), and V g is the speed of sound. The 
irregularity of the winds and the consequent variability of the wind speed may account for 
the fact that, on the occasion of the ASTP experiment, the peak was observed only in two 
revolutions. 

TIDs generated by the supersonic motion of the moon’s shadow across the earth have 
been observed. Some authors have also reported evidence for dawn supersonic terminator 
generation of gravity waves in Langmuir probe data taken by the ESRO 1A satellite. Still, 
few observations have been performed of terminator-related phenomena. BICEPS 
observations may help to bring this subject to the deserved attention of theoreticians and 
experimentalists. 

Topside sounders provide the only means of studying simultaneously details of height 
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and latitudinal variations of the equatorial anomaly. Consequently BICEPS sounder data will 
be used to test the ability of ionospheric models to predict the anomaly behaviour. Recently 
methods have been developed to derive neutral winds from ionosonde measurements of the 
height of the F-layer [Richards, 1991]. While these methods have been applied to 
mid-latitudes where vertical electrodynamic drift can be ignored, their application to low- 
latitude BICEPS data would provide measurements of vertical ionospheric drift and its 
latitudinal variation. 



Figure 1 

At the Day-to-Night transitions are visible very sharp, 
strong, localized enhancements of electron density. This is 
one of the records of electron density measurements obtain- 
ed by the Apollo-Soyuz Test Project (ASTP), with a phase co- 
herent Doppler link established between the two spacecraft. 

The electron density is shown as a function of time, with each 
second corresponding to about 7.8 Km of covered distance. Time 
is given in UT seconds on July 24, 1975. 
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Ill- Broad-band Waveform Degradation due to Ionospheric 
Effects , in Microwave Communications Links 


By moving to higher and higher frequencies, such as to 20 to 30 GHz, it is generally 
believed that a microwave link in transionospheric paths ( typically from geosynchro- 
nous height to ground) is not affected in any appreciable way by the ionosphere. We 
want to show in this Section that this is not the case, and that, by performing 
simultaneous observations with a microwave link from geosynchronous orbit to ground, 
and with an orbiting bistatic sounder such as BICEPS, we could collect new evidence 
on the effect of the ionosphere on the communication performance of the microwave 
link. 

It is known that there are amplitude fluctuation phenomena of a few decibel in links 
at 20 to 30 GHz from geosynchronous orbit to ground, and it would be of great interest 
to perform measurements (by BICEPS) of ionospheric irregularities in the microwave pro- 
pagation path, to explain the observed scintillations. However, even with a totally 
unperturbed ionosphere, just the columnar content of electrons along the extended propaga- 
tion path causes a serious degradation of the radiated waveform, when pulse widths are 
1 to 10 nanosecond ( corresponding respectively to 1 Gigabit/sec and 100 Megabit/sec 
data rates) . These rates are currently either in use or contemplated by advanced commu- 
nications satellites, such as ACTS ( Advanced Communication Technology Satellite) . 

We have performed an analysis to ascertain the extent of the problem. Here are our 
findings. 


The ionosphere is a dispersive medium (velocity of propagation is a function of 
frequency, for e.m. waves) and its phase constant has a very simple 

analytical expression, if we neglect the effect of the geomagnetic field and of colli- 
sions (these simplifying assumptions are acceptable at the frequencies of the ACTS links) 

2* it isJLUr _ UJ 

‘ /£ ^ 

index of refraction = 1/ / - — 

V 

free space wavelength = — - — 


where 


'4 


c = velocity of light 


8 

3 10 m/sec 


Co > = (electron plasma frequency) 


e = electron charge = 1.602 10 -19 Coulomb 
m = electron mass = 9.108 10~ 31 Kg 


N e 


£o= free space permittivity = 8.85 10 -12 F/m 
N = ionospheric electron density, el/m 3 


We have, therefore 


(36*: 
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Equation [1] can be expanded in a Taylor series about the link's angular frequency co 


P>(M — p*o A Qjj> — uj 0 j — B (juj— (Oo) -t- 



Following the treatment by Elliott (1957) , as amended by Knop & Cohn (1962) r 
of pulse waveform degradation in a dispersive medium ( the case considered by these 
Authors was a waveguide) , we can perform in a simple way an initial estimate of the 
pulse degradation to be expected in trans ionospheric paths from geosynchronous 
height to ground ( these are the ACTS paths) . First step is the computation of the 
so-called Elliott parameter a ( a_ = — ^ ^ BL* ) , where T is the length of the 
propagating pulse, B has been defined above, and L is the length of the link ( in 
our case from geosynchronous height to ground) . The second step consists of examin- 
ing the family of curves in Knop & Cohn (1962), that show the various degrees of 
waveform degradation, for various values of the parameter a_ , and of selecting the 
one that matches the value of a_ determined above. A more precise estimate could be 
obtained, but it is not done here, by inserting the value of a. that we have found^, 
into Equation (1) of Knop & Cohn (1962) ,that replaces equation (14) of Elliott (1957) 
that was incorrect. We reproduce in Figure 2 ( and in TABLE I) Knop & Cohn’s family 
of curves. 


For the case of interest to ACTS, we want to ascertain whether or not transionospheric 
propagation from geosynchronous height to ground, at 20 -to-30 GHz, causes an appreci- 
able distortion of the ACTS modulation pulses. We take the width of the pulse T to be 
T = l-to-10 nanosecond, where the value T=10 nanosec is the equivalent pulsewidth 
of the 100 Megabit/sec rate, and T= 1 nanosec applies to the case of the forthcoming 
1 Gigabit/sec rate. We start by computing : s £_ 


where 



BL 



r I \r \ / A/ oUL 

( c f (3| ) J 0 

- = 418.9 sec ^ ( for the frequency of 20 GHz) 


WJJL 


total (columnar) electron content of the entire ionosphere, el/in 
( see Figure 3 ) 


44 




From Figure 3, we find that the largest Total Electron Content for Noon-Time ( at 
the equinoxes, and for a high Mean Sun Spot number) is 10 18 el/m 2 . In the same 
Figure, we read the smallest Content to be 2 10 17 el/m 2 . Consequently, we have : 



« 0.517 10~ 10 

for 

[nM- io 18 

- 0.231 10 -1 ° 

for 

f»dt = 2 10 


el/m 
17 . . 


( maximum) 

( minimum) 


The corresponding values of the parameter * are ; 

4 a t~*. w ^ .02} (max) 




0*009 (min) 


for T = 10 nanosec 


0,2/ (max) 




/ 


for T = 1 nanosec 


0,09 (min) 


We can now enter into the family of distorted curves in Knop s Cohn (1962). We see 
that the case of the 10 nanosecond pulse falls between curve A (the reference curve) 
and curve B . The waveform distortion is limited, but still measurable. The 1 nanosecond 
pulse is, on the contrary, heavily distorted ( we fall between curve C and curve D) , 
and this would create serious problems in achieving the goal of 1 Gigabit/sec data 

rate, if left uncorrected. Correction could be sought through the application of inverse 
filtering. 


We conclude that the effect of the ionosphere on the performance of a microwave 
communications link from geosynchronous orbit to ground is significant and that 
BICEPS measurements would be very useful, if done along the microwave propagation 

path, to establish the correlation between ionospheric parameters and link's 
performance . 
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Figure 


Total electron content of the ionosphere (Noon local time) as a 
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IV- Proposed Experiment with Simultaneous Observations 
with BICEPS and ACTS Geosynchronous Satellite 


A satellite that would be very suitable for the simultaneous measurements 
to be performed by BICEPS, is the ACTS Advanced Communications Technology 
Satellite, could BICEPS be launched in orbit in a few years from today, 
while ACTS is still fully operational. 

ACTS was launched to synchronous height by NASA in September 1993 and is expect- t 

ed to have a lifetime of several years. ACTS is a multibeam, Ka-band communications 
satellite, and is located on the meridian of Denver, Colorado. It uses a 30 GHz 
uplink, and a 20 GHz down-link. Its bandwidth is nominally 1 GHz. The data rate 
is at present 100 Gigabit/sec, while consideration is being given to increase it 
to 1 Gigabit/sec. The multibeam antenna has fixed and hopping spot beams. Several 
ground stations have been already established in Continental US, but the satellite 
has also a steerable antenna that can be aimed, for instance, at Hawaii and at 
Alaska ( see Figures 4 and 5) . It would be advisable to establish a terminal in 
Canada, possibly at very high latitudes, close to the Northern edge of the Earth 
surface that is covered by ACTS. This would be an Experimenter Terminal, that would 
work in association with the BICEPS mission. 


The waveforms that ACTS should transmit for us are : (a) a time-domain short pulse, 
or equivalent waveform , with the shortest possible risetime ( 1 nanosecond would be 
compatible with the stated nominal bandwidth of 1 GHz; however, a few nanosecond wou 
be adequate); (b) a pure CW carrier, 20 GHz to 30 GHz; (c) other waveforms, as required 
to characterize the microwave channel. These waveforms could be transmitted sequentially 
and would measure such fundamental channel parameters as time spread (inclusive of multi 
path), and frequency spread (Doppler shift and spread). Second-order channel functions 
could also be measured, such as channel correlation functions. 


The role of BICEPS would be, fundamentally, to perform measurements that would explain 
in terms of ionospheric phenomena and structures, the propagation- induced non ideali- 
ties, as appearing in the microwave link. The reason is that the propagation paths 
from ACTS to ground are trans-ionospheric and pass through the volume of the ionosp ere 
where BICEPS performs its measurements. For a fraction of many orbits, BICEPS (ei e 
the BOLAS or the WISPRS conf iguration)would perform bistatic ionospheric sounding in the 
volume of interest. In particular, BICEPS would detect and measure E-region and F-region 
electron density irregularities. Travelling Ionospheric Disturbances, large scale an 
small-scale disturbances, spread-F phenomena, terminator-related electron dens y en- 
hancements, aurora phenomena ,and other time-varying, highly dynamical features of the 
ionosphere. These features directly affect the microwave link that passes through, 
are the cause of its time and frequency spreads. BICEPS would give an unprecedented oppor- 
tunity to measure simultaneously cause and effects. With this new y acqu re now e 
on hand, it would be feasible to design a signal processor that could correct for iono- 
spheric effects and make it possible to utilize to its fullest the J ls 

link’s hardware provides ( in the case of ACTS, we could prove that 1 GHz bandwxdth is 
becoming fully usable, with the processor, even at the Northern edge of its coverage). 
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Coverage of ACTS antenna beams and multibeam communications package 
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